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Abstract 

One  of  the  authors  has  studied  the  properties  of  a family  of  Riesz  bases  obtained  by 
perturbing  the  Haar  function  using  B -splines.  Although  these  bases  cannot  be  obtained 
by  multiresolution  analyses,  they  have  other  interesting  properties.  The  present  paper 
discusses  how  a discrete  signal  {ar;Q  < r < TV  ~ 1}  can  be  studied  by  considering 
a suitable  function  of  the  form  /(f)  :=  Ylr=o  arfr{t),  so  that  the  existing  theory  for 
functions  defined  over  a continuous  domain  can  be  applied. 


1 Introduction 

In  what  follows  Z will  denote  the  integers  and  1R  the  real  numbers;  t and  x will  always 
denote  real  variables.  The  support  of  a function  / will  be  denoted  by  supp(/),  its 
quadratic  norm  by  ||/||  and  if  / € L(IR)  its  Fourier  transform  is  defined  by 

/(*):=  [ e~txif{t)  dt. 

</]R  * 

In  [3]  we  found  a family  of  affine  wavelet  Riesz  bases  of  B2(1R),  of  bounded  support  and 
arbitrary  degrees  of  smoothness,  obtained  by  smoothing  the  discontinuities  of  the  Haar 
function  using  B-splines.  Although  these  bases  are  not  orthogonal  they  are  symmetric, 
a feature  that  is  lacking  in  orthogonal  wavelets.  Our  bases  can  be  constructed  so  that 
the  difference  between  the  frame  bounds  (which  are  given  explicitly)  can  be  made  as 
small  as  desired.  In  general,  orthogonal  wavelets  are  represented  by  infinite  series,  and 
for  computational  purposes  values  are  generated  over  a discrete  set  using  the  cascade 
algorithm  [2,  5],  Our  bases,  on  the  other  hand,  are  given  in  closed  form.  We  now  briefly 
describe  how  these  wavelets  are  defined  and  introduce  additional  notation  and  make 
assumptions  that  will  be  used  in  the  subsequent  discussion. 

Let  Nm(t)  denote  the  B-spline  of  order  m (m  > 2)  ([1],  Chapter  4),  X[o,m-i](t)  the 
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characteristic  function  of  [0,  m — 1]), 

m—2  m— 2 

g(t)  :=X[0,m-i](*)  Nm{t-k),  gi{t)  :=g(t-m  + 1),  h(t)  :=  (1/2)  Nm(t-k), 
k=0  k= 0 

and  q(t)  :=  g\(t)  - h(t).  For  0 < 6 < 1/2,  let  a.\  = — a.<z  = — c*3  — ct\  = 2 (m  — 1 )/8, 
Pi  = 2 (m  - 1),  fa  = 2(m  - 1)(1  + <5)/<5,  fa  = -fa  = (m-  l)/8, 

P{i}(t ) ■=  (-l)i~1q{ait+fa),i  = 1,2, 3, 4,  p{5}(t)  :=  -(x[i/2-5  ,1/2)  (i)~X[l/2,l/2+«)W)) 

6 

p{6}(t)  :=  X[o,i/2)W  - X[i/2,i)W,  and  ip(t)  :=  ^pW(f). 

i= 1 

We  will  call  ip  the  perturbed  Haar  wavelet.  In  [3]  we  proved  that  supp  (ip)  C [—5, 1 + <5], 
ip  E Cm_2(l R),  and  that  if  ipj,k(t)  :=  2^2ip(2H  — k),  then  {ipj,k‘,j,k  e Z}  is  a Riesz 
basis,  and  we  provided  explicit  upper  and  lower  frame  bounds.  Moreover,  in  [7]  we 
showed  that  given  a function  /x,  the  wavelet  coefficients  (p,  ipjtk)  can  be  computed  in 
0(N)  steps  (where  N is  the  sample  size),  just  as  in  the  orthogonal  case. 

In  this  paper  we  will  discuss  the  application  of  the  perturbed  Haar  wavelet  to  the 
study  of  discrete  signals.  Let  us  first  look  at  the  orthogonal  case  for  comparison. 

Let  p be  an  orthogonal  wavelet  associated  with  a multiresolution  analysis  {Vj;j  G Z) 
and  a scaling  function  <p,  with  the  caveat  that  the  definition  of  multiresolution  analysis 
that  we  are  adopting  is  that  of  [1]  and  [4],  and  therefore  Vj  C Vj+i,j  € Z,  whether 
other  authors,  like  [2]  and  [5]  assume  that  Vj+X  c Vr  If  a :=  {ar;  0 < r < N - 1}  is  an 
arbitrary  sequence  of  real  or  complex  numbers,  then  this  discrete  signal  is  transformed 
into  a continuous  one  by  considering  the  function  u(t)  :=  ar<P(t  ~ r)- 

The  study  of  the  signal  u(t)  has  two  stages:  the  analysis  stage  consists  in  computing 
the  wavelet  coefficients, whereas  the  synthesis  stage  consists  in  reconstructing  the  signal 
from  the  wavelet  coefficients.  If  Wj  denotes  the  closure  of  the  linear  span  of  the  functions 
Pjtk,  j € Z,  then  the  Wj  are  mutually  orthogonal  and  V0  — ®j<oWj.  Since  v E Vo,  it 
turns  out  that  the  wavelet  coefficients  (i ',Pj,k)  vanish  for  j > 0.  Moreover,  since  v(t) 
has  compact  support,  for  each  j < 0 there  is  only  a finite  number  of  nonzero  wavelet 
coefficients. 

With  the  perturbed  Haar  wavelet  we  face  an  additional  problem:  the  spaces  Wj  are  no 
longer  orthogonal,  and  we  can  therefore  no  longer  assume  that  all  the  wavelet  coefficients 
corresponding  to  positive  values  of  j must  vanish.  Moreover,  we  may  not  even  have  a 
scaling  function:  in  [8]  we  showed  that  if  8 = 2* , where  £ is  a negative  integer,  then  the 
perturbed  Haar  wavelet  ip  that  corresponds  to  this  value  of  <5  cannot  be  generated  by  a 
multiresolution  analysis. 

To  overcome  these  difficulties,  we  proceed  as  follows.  Let  n E Z be  such  that  2"  > 
4(m  — 1),  6{1}(f)  :=X[o,2(m-i))(*M0-  &{2}(0  :=  q(4(m-l)-t),  b(t)  :=  b^(t)  + b^(t), 
fr(t)  arb(2nt  — 4(m  — l)r),  and  f(t)  :=  By  a direct  application  of  [3] 

Lemma  6 we  obtain  the  following 

Lemma  1 The  function  b(t ) has  the  following  properties: 

(a)  supp(b ) C [0,4(m  - 1)],  (b)  b E Cm~2(R),  (c)  b(2(m  - 1))  = 1, 
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(d)  £;b( 0)  = £;b(2(m  - 1))  = i,b(4(m  - 1))  - 0,  1 < k < m - 2, 

(e)  The  total  variation  ofb  does  not  exceed  4 (m  — 1),  (f)  \b(t)\  < 1. 

From  the  preceding  lemma  we  conclude  that  supp(/)  C [0,1],  and  that  the  func- 
tions fr  have  disjoint  supports.  This  implies  that  ||/||2  = ||6||2||a||22-n,  where  ||a||2  := 
\ar\2-  We  will  also  use  the  i\  norm:  ||a||i  :=  l°r|-  Note,  moreover,  that 

/ e Cm-2(1R),  and  that  /(21-”(m  - l)(2r  + 1))  - /r(21-n(m- l)(2r  + l))  = arb{2(m  - 

1))  = Or* 

In  theory,  given  all  its  wavelet  coefficients,  the  function  / can  be  reconstructed  using 
the  frame  algorithm  or  other,  even  faster,  algorithms  [5].  However,  since  there  may  be 
an  infinite  number  of  nonzero  wavelet  coefficients,  the  application  of  such  algorithms 
may  not  always  be  practical.  We  will  adopt  an  approximation  approach.  If  A = A(6,m), 
and  B = B(6,m)  are  respectively  the  lower  and  upper  frame  bounds  of  the  Riesz  basis 
generated  by  ip,  hj ik  :=  { f,ipj,k ),  and  L f :=  then  from  the  error 

estimates  for  the  frame  algorithm  we  know  that  ||L/  - /||  < (( B - A)/(B  + -A))||/||. 
Since,  as  remarked  above,  we  can  make  A and  B as  close  to  1 as  we  want  by  making  8 
sufficiently  small,  we  conclude  that  for  every  e > 0 there  is  a St,  such  that  if  0 < 8 < <5o, 
then  ||L/  - / ||  < e||/||.  To  approximate  / using  the  wavelet  coefficients  it  will  therefore 
suffice  to  approximate  L f by  an  operator  of  the  form 

32 

Ef—  'y 

j=ji  kez 

Observe  that  since  / has  bounded  support,  E f reduces  to  a finite  sum. 

Our  objective  will  be  accomplished  by  showing  that  there  is  a constant  K such  that 

£ Wi*  <K||a||2-hl/2. 

k€  Z 

But  first  we  need  to  prove  five  lemmas,  of  some  independent  interest.  We  begin  with 

Lemma  2 Let  {ai,;k  € Z } and  {bh-,k  € Z } be  increasing  sequences  such  that  at,-  < 
bk- 1 < Ofc+i,  k E Z.  Assume  that  fk  € L2{E)  and  that  supp{fk ) Q [ ak,bk ],  and  let 

; . /:=£*«/*■  7,/ien||/||2<2EfceZIIMI2- 

Proof:  If  r < k - 1 then  bT  < 6fc_2  < whereas  if  r > k + 1 then  ar  > ak+ % > bk. 
This  implies  that  if  r ^ k - 1,  k then  fT(t)  = 0 on  [ofc,  bk],  and  we  readily  see  that 

ll/ll2  <2£  r‘lA(t)|2  = 2£||/fc(0ll2.  □ 

ke  zJak  kez 

Lemma  3 Letu  e L2 (M)  be  a function  with  support  in  an  interval  [a,  6]  with  b-a  < 1 . 
If  j < 0,  then 

X^Ku.^fc)!2  <3IIw|I22;'- 

ke  z 

Proof:  Let  j < 0 be  arbitrary  but  fixed,  and  define  I(k ) :=  supp  {ipj,k)  O [a,b].  Then 
I(k)  C [2 ~j(k  - 8),2~i{k  + 8 + 1)]  D [a,  6].  If  I(k)  = 0 then,  either  2 + 5 + 1)  < a, 
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or  2~3(k  — 8)  > b.  This  implies  that  if  I(k ) 0 0,  then  k £ (23a  — S — 1, 2 3b  + (5).  Since 
the  length  of  this  interval  is  less  than  3,  we  conclude  that  there  are  at  most  three  values 
of  k for  which  I(k)  0 0.  In  other  words,  there  are  at  most  three  values  of  k for  which 
hj.k  0 0.  Since  \if{t)\  < 1,  for  any  such  k we  have: 

|2 


< 


[ u(t)ip(2kt  - k)  dt  < 23  f \u(t)\2  dt  f \ip(2kt  — k)Y 
Ji(k ) y/(fc)  J/(fe) 

(b  — a) 23  [ |u(f)|" 

Jl(k) 


dt 


' dt  <23  llull  . □ 


Lemma  4 Let  a,  /?,  7,  cr  £ 1?,  with  a,  7 0 0,  and  define  c(t ) :=  q(at  + 0),  d(t)  := 
</(-;f  + a),  and 

K = 2 { [25/64  + (25/192)273  (m  - l)4  + (m  - 1)2/1024}  . 

If  j > 0 and  i = 5,6,  then 

(a)  £ |(d)Cj,fc)|2  < 2 (4V^o~2  + 1/3)2 2~J;  (b)  ]T  |<d,P$>|2  ^ (2^  + 1/2)' 2~L 


e-(m-i)«/2  _ ((^sin®^)"*-1 


fcez  fcez 

Proof:  (a)  Prom  [3]  p.  3367  (bearing  in  mind  the  slightly  different  definition  of  the 
Fourier  transform),  we  have 

g(x)  = (i/x)e~(m~1^xi/2 

From  [1]  p.  56  (3.2.16), 

Nm(x)  = e~(l72)m*i[(2/a:)sina:/2]m.  (1.1) 

Let 

s(x)  :=  \{2/x)  sin  x/2]m~1 . 

Then 

gi(x)  = e^m-1)xig(x)  = {i / x)[e-^m~l)xi  - e-(3/2)(m-i)*i s(x)\. 

Since 

1 m— 2 1 — (m— l)a;i 

fr(»)  = 5 E e-^Nm{x)  = (1/2)  1-e_xi  Nm(x), 

k= 0 

a straightforward  computation  yields 

h(x)  = -i/{2x)[e-^2^m-^xi  - e-(3/2)("l-1)^]s(m), 

whence 


q(x) 


This  implies  that 


x 


e ~'~”|cos(m  — l)a;  — s(a;)  cos  -(m  — l)x 


+z(2s(a;)  sin  - (to-  — l)a:  - sin(m  - l)ic)]. 

u 


(^)r  < 8x  , x 0:0. 


(1.2) 

(1.3) 
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On  the  other  hand, 


where 


q(x)~ix~1e  ^ [(vi  + v2)  + i(v3  + n4)] , 


Vx  :=  cos(m  - l)x  — cos(l/2)(m  — l)x,  V2  :=  [1  — s(x)]  cos(l/2)(m  — l)x, 

V3  :=  s(x)  [2  sin(l/2)(m  — \)x  — sin(m  — l)m] , V4  :=  [s(x)  - 1] sin(m  - l)x. 

A McLaurin  expansion  shows  that  |wxj  < (5/8 )(m  — \)2x2.  Since  1 — um~ 1 = (1  — 
u)  YlT=o  uk  and  | sin u | < |u|,  we  infer  that 

|1  - s(x)|  < (m  - 1)|  1 - (2/o;)sina;/2|  = (m  - \){2/x)\x/2  - sinat/2|. 

Since  |u  — sinu|  < |u|3/6,  we  conclude  that  ]1  — s(:r)|  < (m  — l)x2/48.  Thus, 

\vz{x)\  < (m  — l)a:2/48,  and  |z^4 (.-r) | < (m  — l)a;2/48. 

Another  McLaurin  expansion  yields  |t!3|  < (5/24 )(m  — l)3|x|3.  Clearly  |v3|  < 3;  thus 
( V3 1 = |v3(2/,3|w3|1/3  < (25/192)1/3(m  - I)2#2.  Since 

|g(z)|2  = rc_2[(vi  + v 2)2  + (v3  + u4)2]  < 2x~2[v\  + v\  + v\  + u2], 


we  deduce  that 


|g(a:)|2  < Kx2. 


Prom  Plancherel’s  identity  we  have: 


(d,Cj  k)  = 2 j/2  [ d(t)c{2jt  - k)  dt  = 2j'/2/(2tt)  / ekxic{x)d{2H)  dx 

m m 

r2ir  ^ 

= 2^2/(2i r)  / V^c(a:  + 27rr)d(2J(a;  + 27rr))  dx. 

•/o  fc€Z 

This  means  that  {2 ~^2{d,Cjtk)',k  € Z}  is  the  sequence  of  Fourier  coefficients  of  the 
function  J2kez^(x  + 2/Kr)d{2^{x  + 2~r)j.  Thus,  applying  Bessel’s  identity  and  then  the 
Cauchy-Schwarz  inequality  twice  (once  for  sums  and  once  for  integrals),  we  have: 

p2ir  ^ 2 

27t2“:7  K^, cj,fc)|2  = / ^c(x  + 2irr)d(2j(x  + 27rr))  dx 

k€ Z r€Z 

< f J"|c(m)«i(2-2m) | + \c(x  - 27r)d(2;?(x  — 27t))|  + j ^ c(:r  + 27rr)d(2;'(a:  + 27rr))  j dx 

-'°  r#0,-l 

(/■2?r  ^ \ 1/2  / i-2tt  ^ \ 1/2 

J \c{x)d{2^ x)\2  dxj  + yj  \c{x  - 2ix)d{2i{x  - 2-zr))  |2dxj 

0/»27r  ^ \l/2*ii 

| c(x+27rr)d(2i(a:+27rr))|2 da;) 

0 ,*0,-1  ' 

=:  ^\/5i  + %/$2  + \/^3^  • 
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(1.5) 

(1.6) 

(1.7) 


Since  c{x)  = a~1e^^OI'>lctq(a~1x),  (1.3)  implies  that 

\c(x  + 27rr)|2  < 8|x  + 27rr|~2,  x^=2tt r, 

whereas  from  (1.4)  we  see  that 

\c(x  + 2nr)\2  < K a~4\x  + 2irr\2. 

Since  d{x)  = (1.3)  also  implies  that 

|d(2J(x  + 2nr))\2  < 4 ~J+1  2|x  + 27rr|-2,  x ^ 2nr. 

Since  S\  is  obtained  by  integrating  the  product  of  the  left-side  members  of  (1.6)  and 
(1.7)  (with  r = 0)  over  an  interval  of  length  2n,  we  readily  see  that 


A similar  argument  yields 


From  Minkowski’s  inequality 

S3  < ^ |c(x  + 27rr)|2  ^ jd(2J(x  + 27rr)) 

r#  0,-1  r^O.-l 

If  x € [0,  27t]  and  r > 1 then  from  (1.5)  we  have: 

|c(x  + 27rr)|2  < 27r“2  yy  r~2  — 1/3, 


Si  < 1 6nKa~H-j. 

(1.8) 

S2  < 167TAT  a~44~^. 

(1.9) 

dx. 


r>l 


r>  1 


whereas  (1.7)  implies  that 

|d(2J(x  + 27rr))|  < 24“-77r-2  y^r-2  = 4-J /3. 


r>l 


r>l 


Similarly, 


and 


yy  |c(x  + 27rr)|2  < 27r  2y^r  2 = 1/3, 

r<— 2 r>l 

yy  |d(2'7  (x  + 27rr))|  < 2 4~J7r~2  y^  r-2  = 4~J/3, 


r<—2 


r>  1 


whence  we  conclude  that  5,3  < (47r/9)4  j.  Combining  (1.8),  (1.9)  and  the  preceding 
inequality,  the  assertion  follows. 

(b)  Note  that  pi6 1 is  pl5l  with  <5  = 1/2.  Since  pl5l(x)  = 2?x_1e~(1//2)XJ(l  — cosfe), 
we  see  that  ___ . 

|pl5l(x  + 27rr)|2  < 4|x  + 27rr|-2,  x^2'kt.  (1.10) 

On  the  other  hand,  the  inequality  |1  — cos<5x|  < (l/2)< 52x2  implies  that  |pl5I(x)|  < 52\x\\ 

therefore  

|pf5l(x  + 27rr)|2  < 54 \x  + 27rr|2.  (1.11) 
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We  now  repeat  the  argument  employed  in  (a),  using  (1.10)  instead  of  (1.5),  (1.11)  instead 
of  (1.6),  and  bearing  in  mind  that  5 < 1/2.  □ 

We  now  find  bounds  for  the  quadratic  norms  of  q{t)  and  b(t). 

Lemma  5 (a)  ||i/>||  < 1;  (b)  ||6||  < 2(m  - 1). 

Proof:  (a)  [1]  Theorem  4.3  implies  that  the  functions  Nm  are  nonnegative.  This  implies 
that  both  g and  h are  nonnegative.  In  the  proof  of  [3]  Lemma  6(f)  we  show  that 


whence 


f g(t)dt  — f h(t)  dt  — (m  - l)/2, 
J] R VIR 

/ |g(t)|  dt  < m — 1. 

Ji R 


Moreover,  \q(t)\  < 1 ([3]  Lemma  6(h)).  Thus, 


Therefore 


This  implies  that 


[ IqCt)]2  dt  < [ |q(t)|  dt  < m — 1. 

m Ju 

'j 

[ \p^(t)\2dt  = (S/2(m-l))  [ \q(t)\2dt<6/2,  i=  1,2, 3, 4. 
JJR  J1R 

Lies  that 

[ Mt)\2dt  <4<5/2  + f \p^(t)-p^{t)\2dt^2S  + {l-26)  = l. 

Jn  m 

[ \b(t)\2  dt  < [ \b{t)\dt  — 2 f \q(t)\ dt  < 2(m  — 1) . □ 

JM  JlR  J] R 


Theorem  1 

(a)  If  j < 0, 


< 2%/6 (m  - l)||a||  2(j' 


(b)  Let  K be  defined  as  in  Lemma  4-  If  j > 0, 

J2(f,^M,k  < 8 [y/2  (K21~n  + 1/3)  + yfi+  1/3]  ||a||t  2~^2.  . 

kez 

Proof:  Assume  first  that  j < 0.  Applying  Lemma  2,  Lemma  3,  and  Lemma  5,  we  have: 
2 

^2(f>fj,k)i>j,k  <2  \\(f^j,kHj,kf  = \(f^3,k)  i2 

fcGZ  fcEZ  AiEZ 

< 2 £ |(/,Vt*>|2  < 6||/||2  2*  < 6||6||2||a||22^"  < 24(m  - l)2||a||22^". 
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Assume  now  that  j > 0.  Setting  bi^(t)  :=  arb M (2 nt  - 4 (m  - l)r),  we  see  that  fr(t ) = 
bP^ (£)  + bj-2\t).  Thus, 

2 6 N-l  I 


fee  z 


sEEE 

t=l  <=i  r=0 


'j.fe 


fcez 


Applying  Lemma  2 and  Lemma  5 as  above,  we  see  that 


BUMSw*. 


feez 


< 


2E  #}4S> 


feez 


Since  the  Fourier  transforms  of  q(f)  and  X[o,2(m-i))?(£)  are  identical,  and  the  functions 

bP^  are  of  the  form  ar  q(at  + (3)  or  aT  X[o,2(m-i))(at  + P)q(ctt  + p)  with  |a|  = 2",  from 
Lemma  4 we  have: 

E|(,M^r,2K|2(2^  + 1/3)22-,  £ = 1,2,3, 4, 


fee  z 


and 


E |<^}4i}>f  ^ Kl"2  (v^+  1/3)22-j',  * = 5,6, 


feez 


whence  the  assertion  readily  follows. 


□ 
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